Multi-band spectroscopy of inhomogeneous 
Mott- insulator states of ultracold bosons 



D. Clement, N. Fabbri, L. Fallani, C. Fort and M. Inguscio 

LENS, Dipartimento di Fisica, Universita di Firenze and INFM-CNR, via Nello 
Carrara 1, 1-50019 Sesto Fiorentino (FI), Italy. 

E-mail: clement@leiis.uiiifi.it 

Abstract. In this work, we use inelastic scattering of light to study the response of 
inhomogeneous Mott-insulator gases to external excitations. The experimental setup 
and procedure to probe the atomic Mott states are presented in detail. We discuss the 
link between the energy absorbed by the gases and accessible experimental parameters 
as well as the linearity of the response to the scattering of light. We investigate the 
excitations of the system in multiple energy bands and a band-mapping technique 
allows us to identify band and momentum of the excited atoms. In addition the 
momentum distribution in the Mott states which is spread over the entire first Brillouin 
zone enables us to reconstruct the dispersion relation in the high energy bands using 
a single Bragg excitation with a fixed momentum transfer. 
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1. Introduction 



Many-body quantum systems where strong correlations between particles play a central 
role are among the most intriguing physical systems since no simple picture captures 
their behaviour. Examples of this kind can be found in correlated electronic systems such 
as high-Tc superconductors [JJ, quantum liquids such as liquid Helium [2] or interacting 
ID systems [3] as carbon nanotubes [I], quantum wires [5] and organic conductors 
[6]. A large amount of research work has been devoted to the study of such systems 
within different fields of physics, yet many questions have still to be addressed. In 
this prospect new experimental possibilities have been opened by the realization of 
degenerate quantum gases and the recent development of technics to manipulate them 
in optical lattices. Indeed, in the past years ultra-cold atomic gases have been loaded 
in optical lattices to create strongly correlated quantum phases in a highly controllable 
manner. 
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Among these experiments, different atomic insulating phases have been realized 
demonstrating the versatility of gaseous systems. These include bosonic Mott insulators 
in one-dimensional (ID), two-dimensional (2D) and three-dimensional (3D) systems 
[3 EJ [9], disordered bosonic insulating phases [10] as well as fermionic Mott insulators 
[TT1 IT2"] . The characterization of these insulating quantum phases has enlightened 
many of their properties. Experiments have demonstrated the presence of a gap in 
the spectrum of Mott insulators [7J [9] and its vanishing in a disordered insulating 
phase |10j . A noise correlation analysis of time-of- flight pictures has shown the spatial 
order of the atomic distribution [T3] and its controlled modification when modifying 
the lattice potential [15]. Studies of the suppression of compressibility has been 
performed in fermionic [IT], [12] and bosonic [13J Mott insulators. In addition to these 
properties expected in homogeneous insulating phases, the presence of a trap enriches 
the experimental situation while creating alternate regions of Mott with different filling 
factors and superfluid states [TT], [18]. This shell structure has been identified in 
experiments and clearly related to the trapping potential pi| EH]. Recently, inelastic 
light scattering has been used to measure the excitation spectrum of inhomogeneous 
Mott insulating states [21] • This spectroscopic measurement performed at non-zero 
momentum transfer and in the linear response regime gives a direct access to the 
dynamical structure factor S of those complex phases [22j [23], [2U [25] [26] [27] . 

In this paper we report the measurement of the excitation spectrum of ID bosonic 
Mott insulating states over several energy bands. The band structure is induced by the 
presence of the optical lattice which drives the system into a Mott insulating state. A 
particular attention is paid to the details of the spectroscopic technic. We discuss the 
measurement of the energy absorbed by the atomic system when excitations are created 
as well as the linearity of its response. This paper also extends our previous work [21] the 
scope of which was centered on the study of excitations lying in the lowest energy band. 
We present experimental results for the spectra in the higher energy bands showing 
novel experimental signatures related to the one-particle spectral function of the Mott 
state. A band-mapping technic is used to identify the momentum of the excited atoms 
contributing to the different parts of the spectra in the inhomogeneous Mott state. 

The paper is organized as in the following. In section [2] we describe the principles 
and the experimental implementation of inelastic light scattering (Bragg spectroscopy) 
as a probe of excitations of ultra-cold atoms in optical lattices. In section [3] we 
present the experimental sequence allowing to extract information about the amount 
of excitation created by the Bragg beams. We discuss the choice of the experimental 
parameters in order to work in the linear regime, in particular for the response of Mott 
insulating states. Section @] is devoted to the multi-band spectra of the inhomogeneous 
Mott insulating states and to the band-mapping technic. The latter allows us to directly 
construct the dispersion relation in the high energy bands induced by the optical lattice. 
Finally we discuss the relation of the high-energy band spectra with the measurement 
of the one-particle spectral function of the atomic Mott state. 
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2. Inelastic light scattering from correlated gases in optical lattices 

2.1. Inelastic light scattering as a probe of excitations 

The result of weak inelastic scattering of waves or particles by many-body systems may 
be described within the Born approximation and expressed in terms of the dynamic 
structure factor S [28]. The dynamic structure factor is the Fourier transform of 
the density-density correlation function and it carries information on the excitation 
spectrum of the system, independently from the external probe. The knowledge of the 
dynamic structure factor S obtained from inelastic scattering processes has proved to be 
central in describing many-body systems, from the electrons in solids [29] to superfluid 
Helium [2]. 

Inelastic light scattering has been applied to gaseous Bose-Einstein condensates 
to measure the dynamic structure factor [301 [31], [32]. This scattering technic, referred 
to as Bragg spectroscopy, consists in a two-photon transition between two different 
momentum states of the same internal ground-state [30]. The dispersion relation of 
interacting BECs in the mean-field regime [33j [3U [35] , the presence of phase fluctuations 
in elongated BECs [36J as well as signatures of vortices [37] have been investigated using 
this technic. Bragg spectroscopy is also used as a tool to coherently manipulate atomic 
clouds for interferometric schemes [38] or for thermodynamics studies [39, 40J. More 
recently it has succeeded in providing novel information about strongly interacting 3D 
Bose [41] and Fermi [42] gases close to Feshbach resonances as well as correlated ID 
Bose gases across the transition from the superfluid to the Mott insulator state |21j . 

2.2. Two-photon Bragg transitions 

The Bragg spectroscopy technic consists in shining the atoms for a finite time ts with 
two laser beams detuned by 5 from an atomic transition. During the Bragg pulse, atoms 
absorb photons from one beam and are stimulated to emit photons in the second beam. 
The resonance condition for the two-photon transition consists in conserving energy and 
momentum. An atom with an initial momentum p i ends up in the same internal state 
with a final momentum Pf = Pi + hc\ where ftq is the momentum transfer given by the 
two-photon process. The energy difference between the initial and final atomic states 
is also given by the Bragg beams, which are detuned from each other by a frequency 
v — v\ — v-i where V\ and vi are the frequencies of the two beams (see Fig. [1]). In 
practice, by tuning the frequency difference u, one can match the resonance energy for 
the two-photon transition process to occur at a given momentum transfer hq . In other 
words, this corresponds to the condition for which the atomic wave-function is efficiently 
diffracted on the moving optical lattice created by the Bragg beams. We note Vb the 
amplitude of the lattice induced by the Bragg beams and f2s/27r = Vb/2/i the associated 
Rabi frequency. 

The modulus of the momentum transfer ftq given to the atoms is set by the 
wavelength of the Bragg beams and the angle 9 between them, namely qo = 
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47T /\ B sin(0/2) (see Fig.[T]). Note that the direction of q is perpendicular to the bisector 
between the two beams. The detuning v between the two Bragg beams controls not only 
the energy hv transferred to the atoms but also the sign of the momentum transfer. For a 
two-photon transition towards an excited state with higher kinetic energy, atoms absorb 
a photon from the beam with the higher-energy photons and are stimulated to emit a 
photon in the beam with the lower-energy photons. Therefore, with our convention, 
when the detuning v = V\ — v 2 is positive atoms absorb photons from beam 1, emit into 
beam 2 and the momentum transfer is +/iqo (see inset in Fig. [1]). When the detuning v 
is negative atoms absorb from beam 2, emit into beam 1 and the momentum transfer 
is opposite, i.e. —hqo. 




Figure 1. Experimental setup: two laser beams (red), detuned with respect to each 
other by a controllable frequency v, are shone onto the atoms in order to induce a two- 
photon transition. The direction of each beam can be changed thanks to a rotation 
plate and a translation stage that allow a precise control of the angle 8. The axes 
(Ox,Oy) coincide with the axes of the lattice beam (see text below). Inset: parabolic 
energy spectrum of free particles. Depending on the sign of the relative detuning v 
between the Bragg beams atoms absorb photons from beam 1 and emit in beam 2 (for 
v > 0, solid lines) or absorb from beam 2 and emit in beam 1 (for v < 0, dashed lines), 
changing the direction of the momentum transfer hq . 

In our experimental setup, a laser beam coming out from a laser diode at wavelength 
Xb = 780 nm, detuned by 5 ~ 300 GHz from the D2 transition of 87 Rb, is separated 
into two beams, each passing through an Acousto-Optic Modulator (AOM). The two 
AO Ms are locked in phase and allow to control the relative detuning v {y = V\ — u 2 <C 5) 
between the two Bragg beams. After passing through the AOMs, each beam is injected 
into a polarization-maintaining fiber. The out-coupler of each fiber is fixed onto a post- 
holder mounted on a rotation plate and a translation stage close to the cell in which 
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the atoms are manipulated. This setup allows an independent control of the angle of 
each of the Bragg beams with respect to the y-axis (see Fig. [Q. In particular the 
angle 9 between the two Bragg beams can be tuned from 30° to 60°, this range being 
determined by the limited optical access to the glass cell. While A# is precisely known 
in the experiment, the angle 9 cannot be measured geometrically with high precision. 
We use the atoms as a sensor to precisely calibrate q . To this purpose we measure 
the momentum transferred to a 3D BEC (in the absence of any optical lattice) in two 
different ways. 

2. 3. Calibration of the momentum transfer 

In the first calibration measurement, we measure the resonant frequency of the Bragg 
spectrum for the 3D BEC. The resonance frequency depends on the momentum transfer 
hq Q , on the strengh of atom-atom interactions through the dispersion relation and on 
the initial velocity of the BEC center of mass. The dispersion relation E(q) of a trapped 
3D BEC in the mean-field regime can be written as 

E( q ) = ^ DA (q)q 2 +(^) 2 (i) 

where clda(q) corresponds to an effective sound velocity within a Local Density 
Approximation and is related to the mean-field interaction term [51]. The two-photon 
transition induced by the Bragg pulse couples the initial state of the 3D BEC with 
momentum p { = hq i to an excited state with momentum P/ = ft(qj + q )- Thus the 
resonance energy of this process is E(q i + q ) — E(q i ). In the absence of interactions, 
this energy reduces to the usual quadratic dependence of a single-particle spectrum, 
i.e. ft 2 (q.j + q ) 2 /2m — h 2 q^/2m. In order to reduce the effect on evaluating q coming 
from interactions we perform Bragg spectroscopy on dilute 3D BECs after a time-of- 
flight. When the magnetic trap is switched off the BEC acquires a spurious non-zero 
momentum hq { and two unknown parameters have to be determined, namely q and 
qj. We measure the spectra at positive (%q) and negative (—hq ) momentum transfer. 
The two resonance frequencies E(q t + q ) — E(q i ) and E(q i — q ) — E(q i ) allow us to 
precisely determine q in the experiment. We obtain q = 0.97(3)k L where k L = 2n/\ L 
is the wave- vector of the lattice beams (see text below) at the wavelength A^ = 830 nm. 

In the experiment, the bisector of the two Bragg beams is not exactly perpendicular 
to the y-axis and the momentum transfer along the y-axis is the projection Hqo,y — 
4h~7i cos(a) / \b sin(6 l /2) where a is the angle as defined on Fig. [TJ Note that the axes on 
Fig. [T]are defined as the axes of the lattice beams (see below). To measure the angle 
a we diffract the atoms in the Raman-Nath regime [16] with very short light pulses (of 
typical duration 3 [is). In this diffraction regime, many clouds of diffracted atoms are 
observed allowing a good estimate of the axis of the momentum transferred to the atoms. 
We use such a procedure for both the lattice beams and the Bragg beams separately, 
measuring the angle a = 9.5(1)° (see Fig. [Q. As a consequence the projection of the 
momentum transfer along the y-axis is HqQ iV = 0.96(3) hkjj. 
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In a second set of calibration measurements, we use the diffracted atoms of the 3D 
BEC by the moving lattice created with the Bragg beams. By letting the atoms fall 
under gravity after the Bragg pulse for a long enough time-of-flight (£tof = 10 — 30 ms), 
the diffracted atoms separate from the atoms which have not undergone the two-photon 
transition. The distance between the two clouds is frqo.y^TOF /m in the yOz plane where 
absorption images are taken, m being the atomic mass. By fitting the distance between 
the diffracted and non-diffracted atomic clouds as a function of ^tof, we measure 
<io,y = 0.97(4)/cl in good agreement with the previous measurement. 

3. Amount of excitation and linear response regime 

3.1. Energy transfer by the Bragg beams to a gas loaded in optical lattices 

3.1.1. Experimental setup and time sequence We use the Bragg spectroscopic scheme 
described in section [2] to probe the response of correlated ID Bose gases in the linear 
response regime [21] . Correlated Bose phases are created by loading a 3D BEC of 87 Rb 
atoms in a 3D optical lattice as described in previous papers (TO! EI] ■ In brief, a 3D 
BEC of N ~ 1.5 x 10 5 atoms is produced in a magnetic trap which frequencies are 
&x — u z = 2tt x 90 Hz and u y = 2tt x 8.9 Hz. Then it is adiabatically loaded in a 
3D optical lattice created with three pairs of counter-propagating laser beams at the 
wavelength \l = 830 nm. The amplitudes V% of the optical lattice along each axis 
% = x,y,z are expressed in units of the recoil energy E R = h 2 /2m\ 2 L , Vi = SiE R . 
The optical lattices are ramped up to their final values Sj with an exponential ramp of 
duration 140 ms and time constant 30 ms. 

In all the experiments described in this paper, a 2D optical lattice in the plane xOz 
with amplitudes s x = s z = s± = 35 is used to create an array of ID tubes of atoms. At 
this amplitude s± and in the absence of the longitudinal lattice (s y = 0), we measure 
the ratio of the frequencies of the quadrupole mode vq to the dipole mode ud to check 
that the ID gases are Bose-condensed. For holding times in the 2D lattice (at s± = 35) 
up to 50 ms, this ratio is close to the value \/3 corresponding to the expected value 
for ID Bose-condensed gases close to the mean- field regime [4*61 147]. For longer holding 
times, the ratio increases up to the value 2 expected for a thermal ID Bose gas. All the 
experiments are performed with a typical holding time at the amplitude s± = 35 equal 
to 10 ms ensuring the ID gases are Bose-condensed in the absence of the longitudinal 
lattice (s y = 0). 

The amplitude s y of the optical lattice is changed to tune the ID gases from a 
correlated superfluid state (s y = 0) to inhomogeneous Mott insulating states (s y > 6). 
For large amplitudes s y the dynamics of ID gases (when restricted to the lowest energy 
band of the optical lattice along Oy) can be described by the Bose- Hubbard Hamiltonian 



P2], 





<i,j> 
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Here a], a>i are the creation and annihilation operator of one boson at site i and = aja, 
is the particle number operator. The on-site interaction energy is given by U and 
the next-neighbour hopping amplitude by J. £j is the slowly-varying energy offset 
experienced by an atom on site % due to the presence of the harmonic confining potential. 

Once the amplitude of the 3D optical lattice has reached its final values (s± = 35, s y ) 
the Bragg pulse is shone onto the system for a time ts (see Fig. [2(a)). In order to detect 
the amount of excitation induced by the Bragg beams in the correlated gaseous systems 
we follow a procedure similar to that used in [7j and described in [21] that consists in 
measuring a quantity related to the increase of energy in the gas. 

3.1.2. Measurement of the energy transfer After shining the Bragg pulse onto the 
atoms the amplitudes of the 3D lattice are ramped down linearly to s± = s y = 5 in 
15 ms where we let the system rethermalize for 5 ms [see Fig. [2(a)] . Then the magnetic 
trap and the 3D lattice are switched off abruptly and an absorption image of the atomic 
distribution in the yOz plane is taken after a typical time-of-flight £tof — 22 ms. 
Expanding from a phase coherent state in a 3D optical lattice (s± — s y — 5), the atomic 
distribution exhibits an interference pattern which is the analogous of the diffraction 
pattern of light from a grating (see Fig. [2(b) and [13]). From this interference pattern 
we extract the RMS widths o y and a z of the central peak, the increases of which are 
related to that of the energy of the system. Indeed, when the Bragg excitation is tuned 
out of resonance both widths a y and a z are equal to the widths measured in the absence 
of the Bragg pulse while on resonance they increase [see Fig. [2(c)]. Moreover we find 
that o y and a z have the same dependence with the detuning v between the Bragg beams 
as expected from an efficient re-thermalization process in each spatial direction when 
the lattices are ramped down. In the following we will consider an average width a 
defined as o = ^o y o z . 

To be more quantitative on the relation between the increase of cr, noted Act, 
and the energy absorbed by the gas, we make a comparison of Act after a resonant 
Bragg excitation with a measurement of the temperature increase AT realized in the 
same experimental conditions (in particular for an identical Bragg excitation). The 
temperature is obtained by ramping down adiabatically the optical lattices with an 
exponential ramp after the Bragg pulse, letting the system thermalize in the magnetic 
trap for 1 s and measuring the condensate fraction of the 3D atomic cloud [see Fig.[2](d)]. 
The comparison is performed by setting the detuning v on the resonance towards the 
second band in the case of an inhomogeneous Mott insulator state with s y = 13 (see 
Section Hj). By changing the power in the Bragg beams, Act is tuned over a range 
typical of that used in experiments. The results are presented in Fig. [2] where AT is 
extracted from the decrease of the condensate fraction. This measurement confirms that 
Act is proportional to the energy absorbed by the atomic system, i.e. to the number of 
excitations created. 

One may wonder whether the spectra could be obtained measuring the condensate 
fraction instead of the width cr. This could be the case when the amplitude of the 
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Figure 2. (a) Experimental sequence to measure the increase Act of the width of the 
central peak of the interference pattern after a Bragg pulse. The picture (b) [resp. 
(c)] corresponds to an absorption image taken after a non-resonant (resp. resonant) 
Bragg pulse, (d) Experimental sequence to measure the increase of temperature AT 
after a Bragg pulse and a thcrmalization time of 1 second. The picture (e) [resp. (f)] 
corresponds to an absorption image taken after a non-resonant (resp. resonant) Bragg 
pulse, (g) Act as a function of the increase of temperature AT. The points correspond 
to different power of the Bragg beams with a fixed pulse length of 3 ms. 



response of the ID gases is large as in the superfluid regime or in the Mott regime 
for high-energy bands. On the contrary, the amplitude of the response of the Mott 
insulating state within the lowest energy band is low as the system exhibits its insulating 
behaviour. In this case, we observe that the measurement of Act is more sensitive than 
that of the condensate fraction, allowing to detect small excitations. 

3.2. Amount of excitation A 

We now describe the quantity we will be referring to as the amount of excitations. This 
quantity is proportional to the increase Act of the central peak width of the interference 
pattern obtained after a time-of- flight (see Fig. [21(b) and text below). It is rescaled as 
explained in the following. 

The width a of the central peak has a weak dependence on the total atom number N 
which is related to the presence of interactions, independently from the Bragg excitation. 
In the experiment, the total atom number typically fluctuates by 20% from shot-to- 
shot. In order to remove the small variations (up to 8%) on the measurement of a 
coming from fluctuations of N, we calibrate the increase of a with N by changing the 
atom number over a wide range 10 5 — 10 6 . This scaling has been used to substract, 
for each measurement of a, the contribution to a coming from the fluctuations in the 
atom number. To display the spectrum the width in the absence of Bragg excitation 
is substracted in order to plot the increase Act. In addition, we rescale Act by the 
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Figure 3. Bragg spectrum of an inhomogeneous Mott insulating state with lattice 
amplitudes sj_ = 35 and s y = 13. The dashed (resp. dotted) vertical line marks the 
resonance energy A p h ~ U (resp. 2A p h ~ 2U) of the Mott state. The red arrow 
indicates the detuning v at which the linearity of the response has been tested. It 
corresponds to the resonance energy of a particle- hole excitation in the Mott regions. 

parameters of the Bragg excitation, namely by the factor (Pb£b/£) where Pb is the 
power of the Bragg beams, S their detuning with respect to the atomic transition and ts 
the duration of the Bragg pulse [H] . This procedure allows to compare the amplitudes 
of the response for spectra taken with slightly different parameters of the Bragg beams. 
The linear dependence of Act with [Pb^bI S)~ x holds in the linear response regime, the 
validity of which we discuss below. In the following we will refer to the rescaled Act, 
called A, as to the amount of excitation transferred to the system. The plot of A as a 
function of the detuning v between the two Bragg beams gives the excitation spectrum. 
In Fig. [3] we show an example of a spectrum of the inhomogeneous Mott insulating state 
(see text below and |21j). 

3.3. Linear response regime 

The linearity of the response of the system to the Bragg excitation is studied monitoring 
how Act varies with the parameters of the Bragg excitation on resonance. In [2T] we 
have identified the particle-hole excitation energy A p h of the Mott insulating state. Here 
we fix the frequency v = A p ^/h (i.e. on the resonant peak of the Mott insulating state, 
see red arrow on Fig. |HJ) to study the variations of Act as a function of the amplitude 
Vb of the Bragg lattice and the duration At B of the pulse, with the results plotted on 
Fig. HI In the range of parameters used to monitor the spectra (for example Vb = 0.5Pr 
and AtB — 6 ms for the spectrum in Fig. [3] at s y — 13), the increase of Act is linear 
both with Vb and Ate, demonstrating that the energy transferred to the gases by the 
Bragg pulse increases linearly with both those quantities. 
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The inset of Fig.Hb) compares the variation of Aa as a function of the pulse duration 
At B for two amplitudes Vb, namely Vb = 0.5E R (dots) and Vb = 2.0-Er (circles) at 
v = Aph/h for s y = 13. For the larger amplitude V B = 2.0-Er the response to the Bragg 
excitation exhibits a saturation and the time scale of the linear regime is shorter. We 
note that we did not observe Rabi oscillations as in the case of an excitation from the 
ground-state of a 3D BEC in a harmonic trap and a Bogoliubov mode [32] . In the latter 
case the excited state is decoupled from the ground-state, a situation which might not 
be true anymore in the Mott state where strong correlations are present. The increase 
of Aa also saturates for large amplitudes Vb as depicted in the inset of Fig. @Ji). 
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Figure 4. (a) Act as a function of the amplitude Vb of the lattice induced by the 
Bragg beams which is proportional to the power Pb- The measurement is performed 
at a fixed energy transfer v — A p h//i corresponding to the resonance of the particle- 
hole excitation and for a duration of the Bragg pulse Ats = 3 ms. Inset: same 
measurement on larger scale of the amplitude Vb showing saturation for Vb > E R . (b) 
Act as a function of the duration of the Bragg pulse at a given amplitude Vb = 0.5E-R. 
As in (a), the measurement is performed at the energy v — A p h/h of the particle- 
hole excitation. Inset: comparison of the response for two amplitudes Vb, namely 
V B = 0.5E R (dots) and V B = 2.0E R , (circles). 

In addition, we stress that the amplitude Vb of the moving lattice created by 
the Bragg beams for such parameters is much smaller than the amplitude s y of the 
longitudinal lattice: all experiments are performed in a regime where Vb < 0.05V^. 
This regime is different from the one where the Bragg spectrum has been obtained 
using the lattice modulation technic. In the latter case, the amplitude of the lattice is 
typically modulated by ~ 30% due to a momentum transfer q = [3 [TO] ■ In this work, 
the parameters U and J y (see Eq. [2]) describing the system in the Bose-Hubbard model 
are almost unaltered by the additional light potential created by the Bragg beams. 

The probe we use (Bragg beams) weakly affects the state of the system under 
investigation and the amount of excitation varies linearly with the parameters of the 
probe. In this linear response regime the spectra we measure in the lowest energy band 
are proportional to the dynamical structure factor S [28J. 
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4. Excitations in inhomogeneous Mott states towards different energy 
bands 

Bragg spectroscopy has been used in |21j to investigate the response of inhomogeneous 
Mott insulating states within the lowest energy band and to study their properties 
on an energy scale of the order of the Mott gap A p h. Excitations to higher energy 
bands have been recently reported in [18] where the authors focused on the lifetime 
and the coherence of higher bands population. Here, we study the spectrum of 
inhomogeneous Mott insulating states on a large energy scale, focusing on transitions 
to high-energy bands of the optical lattice, for a momentum transferred by the Bragg 
beams q 0)3/ = 0.96(3)^. 

Since the typical energy transferred for excitations in high energy bands is at least 
one order of magnitude larger than the Mott gap (see details later), the excited atoms 
are not correlated to the non-excited ones lying in the lowest energy band. In addition, 
the Bragg excitation being weak, only a small fraction of the atomic cloud is excited. 
Therefore, the small amount of atoms excited in the high energy bands can safely be 
considered as non-interacting particles. This point makes the study of the high-energy 
excitation spectrum of a Mott insulating state much different from that restricted to 
the lowest energy band where atom-atom correlations play a crucial role [2T] . 

4-1. Multi-band spectrum of a Mott insulating state 

As demonstrated for example in [50], the presence of a periodic potential implies the 
existence of different energy bands in the excitation spectrum, the so-called Bloch bands 
[see Fig. 0(a)]. The typical tunneling time of atoms in a lattice is ~ h/J. When it is 
smaller than the time scale of the experiment, a 3D BEC loaded in the optical lattice 
has a long-range phase coherence. Its quasi-momentum distribution in the presence 
of the optical lattice has a small extension around the center of the first Brillouin 
zone qj i2/ = 0. Therefore the resonance in the excitation spectrum is narrow since the 
resonant condition for a two-photon transition between momentum states hqi tV = and 
hqf >y = foqo, y is well defined in energy. Our experimental resolution is good enough to 
observe that this spread in energy is a small fraction of the energy bandwidth. As it is 
shown on Fig. [5(b), the experimental spectrum of this phase-coherent system (s y = 9, 
s± = 0) exhibits several well-defined resonances corresponding to excitations created in 
the different lattice bands [50J. 

The case of the array of ID gases loaded in a ID optical lattice is different, in 
particular when the amplitude of the longitudinal lattice is large enough for the ID gases 
to be in the Mott-insulating regime. On the one hand, strong atom-atom correlations 
are present and are responsible for peculiar features of the Mott state such as the 
existence of an energy gap A p h in the excitation spectrum. On the other hand, the 
quasi-momentum distribution of a Mott insulator is spread over the first Brillouin zone 
(—k L < q i y < k L ) as a consequence of the spatial localization of atoms |52j. Deep 
enough in the Mott regime, this quasi-momentum distribution is expected to be close to 
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homogeneous over the interval [—&£,; ki]. Therefore two-photon transitions towards the 
excited bands are possible from any initial quasi-momentum and Bragg spectroscopy 
can create excitations on a large energy interval, of the order of the energy bandwidth 
of single-particle. In this picture, fixing the frequency between the Bragg beams allows 
one to selectively excite (and therefore address) only a fraction of the atomic cloud 
through the resonance condition: for a given energy transfer lying in the energy interval 
of an excited band, one can always find a populated quasi-momentum state matching 
the resonance condition for being excited. 
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Figure 5. (a) Energy bands of single particles in a periodic potential of amplitude 
s = 9. The grey areas cover the entire energy distribution of each band and they 
are reported in Figs. EJb)-(c). (b) Bragg spectrum of a 3D BEC loaded in a ID 
optical lattice at the amplitude s y — 9 (blue dots) and for a momentum transfer 
Sqo,j, = 0.96(3)?ifcL- The solid blue line is a fit with three Gaussian functions, (c) 
Bragg spectra of an array of ID BECs (sj_ = 35) loaded in a ID optical lattice of 
amplitude s y = 9 (red dots), i.e. the ID gases being in the Mott-insulating state, and 
for a momentum transfer ftqo,j, = 0.96(3)7jfc£. The red solid line is a guide to the eye. 



Fig. [5](c) depicts a spectrum of ID gases in the Mott insulating state (s± = 35, s y = 
9). The energy scale is identical to that of a 3D BEC loaded in a ID optical lattice along 
the y-axis with the same amplitude s y = 9 [Fig. EJ^b)]. On this energy scale, the spectra 
of ID Mott insulating states exhibit several large resonances that have to be identified 
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with transitions towards different energy bands of the optical lattice. The experimental 
measurements are reported with single-particle energy bands as grey areas. The lowest- 
energy band [y < 10 kHz) is much larger than the single-particle energy band due to the 
presence of atom-atom correlations. A detailed analysis of this effect has been the object 
of [21]. The energy transfer corresponding to transitions towards high energy bands is 
tens of kHz much larger than the Mott gap A p h ~ h x 2 kHz [21] . Identifying the final 
momentum state of atoms giving contribution for the different parts of the excitation 
spectrum of Fig. [5b) should allow us to determine to which band the excitations belong 
to. To this purpose we have used a band mapping technic [13] that we describe in the 
following paragraph. 

4-2. Measurement of the momentum of the excited atoms by the Bragg beams 
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Figure 6. (a) Experimental time sequence of the band-mapping technic, (b) Band 
population in the case of a 3D BEC in the presence of a ID optical lattice of amplitude 
s y = 9 after a Bragg pulse with detuning v = 38 kHz. (c) Band population in the case 
of an array of ID BECs loaded in a 3D optical lattices of amplitudes [s± = 35, s y = 9) 
after a Bragg pulse with detuning v = 34 kHz. The horizontal scale is normalized to 
the momentum hk^ of the longitudinal optical lattice in both pictures. 

4-2.1. Band population in a ID optical lattice In order to identify the band towards 
which atoms are excited by the Bragg pulse we switch off abruptly the magnetic trap 
and we ramp down the lattice on a time scale (~ 2 ms) too short for interaction-induced 
thermalization but adiabatic with respect to the atomic motion in a single lattice well 
[4~3] . Doing this, one ensures to preserve the band population. After a time-of-flight, the 
atomic distribution reflects the band population in the momentum space. In practice, 
we ramp down the optical lattice with an exponential ramp of duration 2 ms and time 
constant 0.5 ms as sketched in Fig. [6](a). We first apply this technic to the 3D BEC 
loaded in a ID optical lattice with amplitude s y = 9 corresponding to the spectrum of 
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Fig. [5(b). The density profile along the y-axis after time-of- flight is shown on Fig. E](b) 
when a resonant Bragg pulse creates excitations in the third energy band [y = 38 kHz). 
The central peak around q = corresponds to the non-excited cloud and it is the 
only feature left when the Bragg pulse is tuned out of resonance. On resonance a 
small lateral peak appears corresponding to excited atoms [53] . The momentum of the 
excited atoms lies between 2hk L and 3hkL, proving that these atoms are excited in 
the third band of the optical lattice. Fitting the position of the diffracted cloud with 
respect to the non-diffracted one with a Gaussian function we have measured momentum 
transfers of and 2.8(l)hk L for a Bragg pulse corresponding respectively to the 

transition towards the second and third energy band of the lattice. These results are in 
good agreement with the expected values l.QA{3)hk L and 2.96(3)hk L corresponding to a 
momentum transfer hqo, y = 0.96(3)^^ given by the Bragg beams in our configuration 

4-2.2. Band population in the Mott insulating states We turn to the case of the ID 
gases in the inhomogeneous MI state (s± = 35, s y — 9). The same experimental technic 
as described above is used to map the band population also in this case. The density 
profile obtained from the band mapping of the array of ID gases is shown on Fig. [6](c) . 
The detuning between the Bragg beams is v = 34 kHz and corresponds to an excitation 
that lies in the wide and flat region of the spectrum. As in the case of Fig. [6(b), 
the diffracted atoms are clearly visible and lie in the third band of the optical lattice 
(2/c L < q < 3k L ). We note that the central (non-diffracted) peak is larger and flatter 
over the lowest lattice band with respect to the case of Fig. [6(b). This observation 
indicates that the quasi-momentum distribution of the inhomogeneous MI state covers 
almost uniformly the first band as expected for a system that is not phase- coherent. 

By applying the same fitting procedure as in the case of Figl6(b), we extract the 
momentum hq to the diffracted atoms in the case of a 3D lattice with amplitudes 

= 35, s y ). We repeat this procedure over the entire frequency range corresponding 
to the transition towards the third and fourth bands, i.e. varying the relative detuning v 
of the Bragg beams from 27 to 45 kHz, with the results plotted in Fig.[7]for s y = 10. The 
experimental resolution ~ 400 Hz of the Bragg spectroscopy allows us to do a precise 
mapping on this large energy scale. In Fig. [7] we plot the energy transfer hv given by 
the Bragg beams as a function of the momentum of the excited atoms measured using 
the band-mapping technic. These results are compared with the dispersion relation 
of single-particles in the presence of a periodic potential with an amplitude s y = 10 
showing a good agreement [57J. This demonstrates that: (i) the excitations observed 
over a large energy scale between 27 and 36 kHz (resp. between 37 kHz and 45 kHz) 
correspond to transition towards the third (resp. fourth) energy band of the optical 
lattice ; (ii), the quasi-momentum distribution of the inhomogeneous MI state extends 
over the entire lowest lattice band since the entire energy band can be mapped using a 
Bragg excitation with a fixed momentum transfer frqo, y . 
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Figure 7. Energy dispersion as a function of the quasi-momentum from the band- 
mapping technic used after a Bragg pulse. The dots correspond to the experimental 
measurements and the red and blue colors are associated respectively to excitations in 
the third and fourth band of the optical lattice. The solid lines are the band dispersion 
relation for single particles at s y = 10 and the shaded areas corresponds to the band 
dispersion relation taking into account the experimental uncertainty on the calibration 
of s y (10%). 

4-3. High-energy bands, towards novel information about the Mott state 

Using a band-mapping technic we have demonstrated that the different parts of the 
spectrum measured in a Mott insulating state [see Fig. Etc)] can be attributed to 
the different energy bands induced by the optical lattice. At high-energies, e.g. for 
excitations into the third and the fourth bands, the frequency range of the Mott 
excitations corresponds to the bandwidth of the single-particle spectrum, a feature 
which is expected since in that situation the typical atom-atom correlation energy A p h 
is much smaller than that of the excited atoms. Yet, as we discuss it below, the peculiar 
lineshape of the response within a given band has to be related to the properties of the 
many-body state which is probed. Therefore the identification of the energy scale of 
a given excited band opens the possibility to study the modulation of the response to 
the Bragg pulse within a single band. For example, the response of the Mott insulating 
state to the Bragg spectroscopy shows a structure in amplitude within the third band 
(from 25 to 36 kHz) on Fig. [5tc). In particular, we observe an increased response at 
both edges of the band. A peak is mostly evident close to the higher-energy edge, the 
position of which corresponds to a Bragg transition from an initial quasi-momentum 
qi = 0. We now schematically discuss a way to interpret these features. 

We expect the atoms excited by the Bragg beams in high-energy bands not 
to interact with the non-excited ones since their typical energy (corresponding to 
a frequency of several tens of kHz) is much larger than the particle-hole excitation 
energy A p h ~ h x 2 kHz. In addition, as already mentioned, with our experimental 
parameters only a minor fraction of the atoms is excited. Therefore the transition 
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induced by the Bragg beams in high-energy bands can be considered to happen between 
an initially strongly correlated state (Mott state) and single-particle states: the two- 
photon transition creates a hole in the Mott state and it populates highly energetic 
single-particle states. In other words, the response of the gas to the Bragg pulse 
gives access to the one-particle spectral function of the Mott and it implies both the 
dispersion relation of the hole in the Mott phase and the density of states (DOS) to 
populate the excited band from the initial correlated phase [58]. Information about 
the one-particle spectral function of solid-state systems has proved to be crucial to 
understand properties of correlated states such as high-Tc superconductors [59] and 
recent proposals have demonstrated the interest of measuring the one-particle spectral 
function for ultracold fermions [60]. We anticipate that the information about the one- 
particle spectral function obtained from the high-energy bands spectra might shed new 
light on the bosonic Mott state, in particular close to the SF-MI transition. To our 
knowledge, no theoretical predictions of such an experimental signal measured with 
ultracold bosons exist. 

The experimental situation is more complex than the case of a homogeneous 
Mott state. Due to the presence of the longitudinal trap, the ID gases are in 
the inhomogeneous MI state formed of Mott regions with different fillings separated 
by superfluid areas [18]. This state can be viewed as constituted of domains with 
narrow quasi-momentum distributions (superfluid regions) and domains with almost 
homogeneous quasi-momentum distributions (MI regions), resulting in a non-trivial 
momentum distribution of the global state. As U/ J increases and the system is driven 
deeper in the Mott phase, the number of atoms which do not belong to a Mott state 
drops. For this reason we expect the contribution of the superfluid domains in the 
response to the Bragg excitation not to be significant. However, a complete picture 
taking into account contributions from the different Mott and superfluid regions has to 
be developed to interpret the details of the spectrum of Fig. [5b). We hope that the 
measurements presented in this paper could stimulate such theoretical work. 

5. Conclusion 

In conclusion, we have measured the response of bosonic Mott insulating states to 
excitations induced by inelastic scattering of light. As already demonstrated in [2T] . this 
spectroscopic technic is a promising tool to characterize strongly interacting quantum 
phases. We have presented here in detail the experimental setup and the way the 
measurement has been carried out. In particular we have demonstrated that the quantity 
we measure is proportional to the increase of energy in the gaseous system and that the 
response to the light scattering is in a perturbative and linear regime. 

We have extended our previous work [21] to the study of excitations towards several 
energy bands in Mott insulating states. In contrast to the case of a 3D BEC loaded 
in a ID optical lattice, the spectra of these Mott states exhibit broad resonances in 
energy. The use of a band-mapping technic after applying the Bragg pulse allows us to 
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identify these resonances with transitions towards the different energy bands induced 
by the optical lattice. We also give direct experimental evidence that the momentum 
distribution of Mott states spreads over the entire first Brillouin zone. This property 
enables us to reconstruct the dispersion relation of the high energy bands using a Bragg 
excitation at a fixed momentum transfer. 

Finally the amplitude of the response of the inhomogeneous Mott insulating state 
in high energy bands exhibits peculiar structures within a single excited band. We 
expect that these could give important information on the experimental system and we 
think that this point should be the object of future investigation, both experimental 
and theoretical. 
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